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ABSTRACT 

We report Fermi - LAT observations of the radio-loud AGN SBS 0846+513 (z=0.5835), 
optically classified as a Narrow-Line Seyfert 1 galaxy, together with new and archival 
radio-to-X-ray data. The source was not active at y-ray energies during the first two 
years of Fermi operation. A significant increase in activity was observed during 2010 
October-2011 August. In particular a strong y-ray flare was observed in 2011 June 
reaching an isotropic y-ray luminosity (0.1-300 GeV) of 1.0 xlO 48 erg s _1 , comparable 
to that of the brightest flat spectrum radio quasars, and showing spectral evolution in 
y rays. An apparent superluminal velocity of (8.2±1.5)c in the jet was inferred from 
2011-2012 VLBA images, suggesting the presence of a highly relativistic jet. 

Both the power released by this object during the flaring activity and the apparent 
superluminal velocity are strong indications of the presence of a relativistic jet as pow- 
erful as those of blazars. In addition, variability and spectral properties in radio and 
y-ray bands indicate blazar-like behaviour, suggesting that, except for some distinct 
optical characteristics, SBS 0846+513 could be considered as a young blazar at the 
low end of the blazar’s black hole mass distribution. 

Key words: galaxies: active - galaxies: nuclei galaxies: Seyfert - galaxies: individual 
(SBS 0846+513) - gamma rays 


1 INTRODUCTION 

Active Galactic Nuclei (AGN) are the most luminous persis- 
tent sources of high-energy radiation in the Universe. How- 
ever, only a small percentage of AGN are radio-loud, and 
this characteristic is commonly ascribed to the presence of 
relativistic jets, roughly perpendicular to the accretion disk. 
Accretion of gas on to the supermassive black hole (SMBH) 
is thought to power these collimated jets, even if the nature 
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of the coupling between the accretion disc and the jet is 
still among the o utstanding open q uestions in high-energy 
astrophysics le.g.. iBlandford IkOPOnMeier IllOO.'jl ). Certainly 
relativistic jets are the most extreme example of the power 
that can be generated by a SMBH in the centre of an AGN, 
with app arent bolometric luminosities up to 1 Q 49 ~ 5 0 erg 
s -1 fe.g.. lAckermann et al.|[20ld : iBonnoli et al.l 1201 il l . and 
a large fraction of the power emitted in y rays. 

Before the launch of the Fermi satellite only two 
classes of AGN were known to generate these structures 
and thus to emit up to the y-ray energy range: blazars 
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and radio galaxies, both hosted in giant elliptical galaxies 
(IBlandford fc Rees I11978T) . The first two years of observa- 
tions by the Large Area Telescope (LAT) on board Fermi 
confirmed that these two are the most numero us classes o f 
identified sources in the extragalactic 7-ray sky (lAbdo et al.l 
l2010al : lNolan et a! .1120121 ). but the discovery of variable 7-ray 
emission from 4 radio- loud Narrow-Line Seyfert 1 galaxies 
(NLSls) revealed the presence of a possible emer ging third 
class of AGN with relativistic jets (lAbdo et al.ll2009aj lbl!ch. In 
addition, Fermi - LAT observations gave us the opportunity 
to study in more detail particular subclasses of the estab- 
lished ty pes of 7-ray emitting AGN (e.g. Broad L ine Radio 
Galaxies, lAbdo et al.ll2010bl : iKataoka et al.ll201ll h On the 
contrary, no ra dio-quiet Seyfert galaxie s were detected in 7 
rays until now dAckermann et aljl2012lb 


NL S1 are a cl ass of AGN identified 

by lOsterbrock fe Pogge I (|l985l 4 and characterized by 
the following optical properties: narrow permitted lines 
(FWHM (H,0) < 2000 km s -1 ) emitted from the broad line 
regio n, [OIIIJ/H/3 < 3 (a criterion introduced by Goodrichl 
1 1989h . and a bump due to Fe II (see e.g. IPogge I l2000l . 
for a review). They also exhibit strong X-ra y variability, 
steep X-ray spectra (ph oton indices Fx > 2; iBoller et al.l 
199(|;; Gru]oeet_al. |2010f ). and substantial soft X-ray excess 

I Boiler et al.lll99£ ). These characteristics point to systems 
with smaller mas ses of the central black hole (10 6 -10 8 Mq, 
lYuan et alj 120081 ) than in blazars and radio galaxies and 
higher accretion rates (close to or above the Eddington 
limit, lYuan et~aH l200ij ). NLS1 are generally radio-quiet 

( R < 10, radio-loudness R being defined as ratio of rest- 
frame 1.4 GHz and 4400 A flux densities), w ith only a small 
fract ion classified as radio-loud (< 7%, iKomossa et al.l 
l2006lL even more sparse (~2.5%) are very radio-loud 
NLSls (R > 100), while generally 10%-15% of quasars 
are radio-loud and very radio-loud. In the past, several 
authors investigated the peculiarities of radio- loud NLSls 
with non-simultaneous radio to X-ray data, suggesting 
similarities with t he young stage of quasars or di f ferent 
types of blazars ( Komossa et al.l 120061 : lYuan et ~ahl l200Sl : 
iFoschini et al.ll2009l ). The strong and variable radio emis- 
sion, and the flat radio spectrum together with variability 
studies suggested the presence in some radio- loud NLSls 
of a relativistic jet, confirmed in some objects by the 
Fermi - LAT detection in 7 rays dAbdo et al.l 2009c]). This 
finding poses intriguing questions about the nature of these 
objects, the onset of production of relativistic jets, and the 
cosmological evolution of radio-loud AGN. The impact of 
peculiar characteristics of the central engines in radio-loud 
NLSls, which seem quite different of those of blazars and 
manifest in their peculiar optical characteristics, on the 
7-ray production mechanisms is currently under debate. 

I11 June 2011 high 7-ray activity from SBS 0846+513 
was obse rved by Ferrm-LAT. Pre liminary results were re- 
ported in lDonato fe Perkinsl d201ll 'l. That was an important 
confirma tion of the detec tion of a new 7-ray NLS1 after the 
claim by IFoschini I (1201 ll fl. 


1 The analysis presented in IFoschini I d201 il l was performed with 
the P6_V3 IRFs, starting from the first Fermi-LAT catalogue 
(1FGL) as a reference for the background sources, and over the 
period August 2008 -February 2011. 


First identified with a BL Lac object at redshift z = 1.86 
during a bright state bv lArp et al.l d 19791 '). SBS 0846+513 is 
a high redshift object positioned very clos e to a l ow r edshift 
galaxy at the end of a chain of five galaxies. I Arp et all (Il979l) 
noted also a small nebulosity close to SBS 0846+513 that 
could be a normal red galaxy observed by chance near the 
object due to projection effects. The Sloan Digit al S ky Sur- 
vey (S PSS) spectrum of the source reported in IZhou et all 
d2l)lirj) is typical of NLS1: FWHM (H/3) = (1710 ± 184) 
km s _1 , [OIII] /H/3 ~ 0.32, a nd a strong Fe II bump. Sim- 
ilar values were obtained by lYuan et al.l d2008l ) analysing 
the same SDSS spectrum: FWHM (H/3) = (1810 ± 191) km 
s' * 1 * 11 and [OIII] /H/3 ~ 0.31. Moreover, the SDSS spectrum 
showed that its true redshift is z = 0.5835. iNottalel (ll986l) 
and others suggested that it was a gravitationally lensed 
quasar with variability due to gravitational amplification by 
a star in an intervening galaxy. However, with broader wave- 
length coverage and high resolution, no signs of any inter- 
vening galaxies have been found in the SDSS spectrum. The 
Hubble Space Telescope (HST) image, collected when SBS 
0846+513 was in a faint state (V ~19.7 mag) does not show 
signific ant resolved struct ure, thus no indication of the host 
galaxy dMaoz et al.lll993h . Hig h polarization (> 10%) was 
det ected in optical and radio bv lMoore fe Stockman I (Il98ll) 
and lSitko et al.l 1 19841) . A remarkable optical variability was 
observed in the past with an amplitude of A Y ~ 5 m a g ove r 
1 year and AY ~ 4 mag over 1 month (|Arp et al.l ll979l l . 
The high pola rization , ve ry hig h brightness temperature 
(' T b > 10 13 K: IZhou et ZII200 3), and the large-amplitude 
variability observed in optical were clues for the presence of 
an at least mildly relativistic jet, now confirmed by the LAT 
detection in 7 rays and superluminal motion revealed by our 
new radio observations. 

In this paper we present the detection of SBS 0846+513 
in 7 rays and discuss its characteristics by means of the avail- 
able radio-to-7-ray data. Preliminary results about the high 
re soluti o n VLBA obse rvations of the source were presented 
in lOrienti et al.l (120121 ). The paper is organized as follows: 
in Section 2 we report the LAT data analysis and results. 
In Section 3 we report the result of the Swift observations 
performed in August-September 2011. Radio data collected 
by Medicina, OVRO, Effelsberg, VLBA, and VLA are pre- 
sented in Section 4 and discussed in Section 5 and 6. In 
Section 7 we present the spectral energy distribution (SED) 
modeling, while discussion and concluding remarks are pre- 
sented in Section 8. 

Throughout the paper, a A-CDM cosmology with H 0 
= 71 km s _1 Mpc -1 , 9 a = 0.73, and 9 m = 0.27 is adopted. 
The corresponding luminosity distance at a = 0.5835 is d l = 
3406 Mpc, and 1 arcsec corresponds to a projected distance 
of 6.584 kpc. 


2 FERMI - LAT DATA: SELECTION AND 
ANALYSIS 

The Fermi - LAT is a 7-ray telescope operating from 20 MeV 
to > 300 GeV. The instrument is an array of 4 x 4 identical 
towers, each one consisting of a tracker (where the photons 
are pair-converted) and a calorimeter (where the energies of 
the pair-converted photons are measured). The entire instru- 
ment is covered with an anticoincidence detector to reject 
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Figure 1 . Monthly integrated flux (E > 100 MeV) light curve of 
SBS 0846+513 obtained from 2010 August 4 to 2011 August 4. 
Arrows refer to 2-cr upper limits on the source flux. Upper limits 
are computed when TS < 10. The horizontal line indicates the 
period considered for the SED. 


the charged-particle background. The LAT has a large peak 
effective area (~ 8000cm 2 * for 1 GeV photons), an energy 
resolution typically ~10%, and a field of view of about 2.4 
sr with an angular resolution (68% containment angle) bet- 
ter than 1° for energi es above 1 GeV . Furt her details about 
the LAT are given by I Atwood et all (l2009h . 

The LAT data reported in this paper were collected over 
the first 40 months of Fermi operation, from 2008 August 4 
(MJD 54682) to 2011 December 4 (MJD 55899). During this 
time the LAT instrument operated almost entirely in survey 
mode. The analysis was performed with the ScienceTools 
software package version v9r23pl. The LAT data were ex- 
tracted within a 15° Region of Interest (Rol) centred at 
the radio location of SBS 0846+513. Only events belong- 
ing to the “Source” class were used. In addition, a cut on 
the zenith angle (< 100°) was also applied to reduce con- 
tamination from the Earth limb 7 rays, which are produced 
by cosmic rays interacting with the upper atmosphere. The 
spectral analysis (from which we derived spectral fits and 
photon fluxes) were performed with the post-launch instru- 
ment response functions (IRFs) P7S0URCE_V6 using an un- 
binned maximum likelihood method implemented in the Sci- 
ence tool gtlike. 

The background model used to extract the 7- 
ray signal includes a Galactic diffuse emission compo- 
nent and an isotropic component. The model that we 
adopted for the Galactic component is given by the file 
gal_2yearp7v6_v0.fits, and the isotropic component, which is 
the sum of the extragalactic diffuse emission and the resid- 
ual charged particle background, is parametrized by the file 
iso_p7v6source.tx10. The normalizations of both components 


in the background model were allowed to vary freely during 
the spectral point fitting. 

We examine the significance of the 7-ray signal from 
the sources by means of the Test Statistics (TS) based 
on the likelihood ratio test. The Test Statistic TS = 
2Alog(likelihood) between models with and without the 
source is a measure of the probability of having a 7-ray 
source at the localization specified, which compares models 
whose parameters have been adjus ted to maximise the likeli- 
hood of the data given the model dMattox et al.lll996l ). The 
source model used in gtlike includes all the point sources 
from the second Fermi - LAT catalogue (2FGL; iNolan et al.l 
l2012ll that fall within 20° from the source. The spectra of 
these sources were parametrized by power-law functions, ex- 
cept for 2FGL J0920. 9+4441, and 2FGL J0957.7+5522, for 
which we used a log-parabola in their spectral modeling as 
in the 2FGL catalogue. We removed from the model the 
sources having test statistic TS < 25 and/or fluxes below 
1. Ox IIP 8 photons cm” 1 s” 1 over 40 months and repeated 
the fit. Thus a final fitting procedure has been performed 
with the sources within 10° from SBS 0846+513 included 
with the normalization factors and the photon indices left 
as free parameters. For the sources located between 10° and 
15° we kept the normalization and the photon index fixed 
to the values obtained in the previous fitting procedure. The 
Rol model includes also sources falling between 15° and 20° 
from the target source, which can contribute to the total 
counts observed in the Rol due to the energy depended size 
of the point spread function of the instrument. For these ad- 
ditional sources, normalizations and indices were also fixed 
to the values of the 2FGL catalog. 

SBS 0846+513 was not in the 1FGL or 2FGL cata- 
logues, indicating that the source was not detected with TS 
> 25 in either one year or two years of Fermi observations 
llAbdo et al.ll2010al ; INolan et al. 120121 1. Integrating over the 
first two years of Fermi operation the fit yielded a TS = 
14. The 2-cr upper limit over the first 2 years of LAT data 
is 8.5xl0” 9 photons cm -2 s” 1 in the 0.1-300 GeV energy 
rang<0, assuming the photon index is T = 2.3, which is the 
average value found for low-synchrotron- peake d blazars in 
the Second LAT AGN Catalog (2LAC) llAckermann et ah] 
2012). On the contrary, the fit with a power-law model to 
the data integrated over the third year of Fermi operation 
(2010 August 4 - 2011 August 4; MJD 55412-55777) in the 
0.1-300 GeV energy range results in a TS = 653, with an 
integrated average flux of (6.7 ± 0.5) xlO” 8 photons cm -2 
s” 1 and a photon index F = 2.23 ± 0.05. 

Finally, over the period 2011 August 4-December 4 
(MJD 55777-55899) only a TS = 8 was obtained, with a 2-cr 
upper limit of 1.9xl0” 8 photons cm” 2 s” 1 in the 0.1-300 
GeV energy range, assuming the photon index is F = 2.23. 
During August-December 2011 the source was not detected 
with TS > 10 using a time bin of 1 month, indicating a 
significant variability in the 7-ray activity of the source on 
month time-scales. For the rest of the paper we focused on 
the 7-ray data collected over 2010 August 4-2011 August 4. 

The 7-ray point source localization by means of the 


2 http:/ /fermi. gsfc.nasa.gov/ssc/data/access/lat/Background 

Models.html 


3 To obtain a convergence we performed the fit considering only 
the sources within 10° from SBS 0846+513 and E > 300 MeV 
and extrapolated the flux down to 100 MeV. 
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gtfindsrc tool applied to the 7 rays extracted during the 
third year of observation results in R.A. = 132.48 deg, Dec. 
= +51.13 deg (J2000), with a 95% error circle radius of 
0.06°, at an angular separation of 0 . 01 ° from the radio posi- 
tion of SBS 0846+513 (R.A. = 132.492 deg, Dec. = +51.141 
deg, J2000). This implies a strict spatial association with 
the radio coordinates of the NLS1 SBS 0846+513. The pro- 
posed association could be questioned to some extent due to 
the presence of a low redshift star-forming galaxy almost on 
the line of sight toward SBS 0846+513. And indeed, several 
members of the Local Group and also some nearby starburst 
galaxies have been recently detected by Fermi - LAT and thus 
established as sources of high energy photons. Their 7 -ray 
luminosities range however between 10 40 and ~3 xlO 43 erg 
s -1 , and no variability in th eir 7 -ra y emission was neither 
observed nor expected (lAbdo et alj|2010eh . The very large 
7 -ray luminosity of the LAT object discussed here and its 
flaring activity therefore rule out the foreground galaxy as 
the source of the detected 7 rays. 

In order to test for curvature in the 7 -ray spectrum 
of SBS 0846+513 two alternative spectral models with re- 
spect to the power-law (PL) were used: the log parabola 
(LP), dN/dE oc E / E -°‘-P l °z( E / E b) dLandau et al.1 1 19861: 
iMassaro et al.1 [20041. and the broken power law (BPL). In 
the case of the log parabola the parameter a is the spectral 
slope at the energy Eq and the parameter /? measures the 
curvature around the peak. We fixed the reference energy 
E 0 to 300 MeV. We used a likelihood ratio test to check the 
PL model (null hypothesis) against the LP model (alterna- 
tive h ypothe sis). T hese values may be compared, following 
I Nolan et al.l (2012t h by defining the curvature Test Statis- 
tic TS CU rve=2(log LLP-log Lpl)=10 corresponding to ~3.3-<r 
difference. This value is below the threshold of TS CU rve=T 6 
applied in lNolan et al.1 (l2012t l for defining a significant cur- 
vature, thus suggesting only a hint of spectral curvature in 
the average 7 -ray spectrum of SBS 0846+513. We tested also 
the broken power-law model, with photon indices T 1 below 
and P 2 above the break energy Inbreak- We fixed inbreak at 1.4 
GeV. This value was estimated studying the profile of the 
Likelihood function, fixing inbreak at different values between 
100 MeV and 3 GeV with a step of 50 MeV and finding the 
minimum log(likelihood) value. Also for this spectral model 
we obtain TS CU rve=2(log LspL-log Lpl)=10, therefore we 
adopt the power-law for the following LAT analysis. The fit 
results are reported in Table [TJ 

Figure [T] shows the 7 -ray light curve of the third year 
of Fermi observations built using 1-month time bins. For 
each time bin the photon index of SBS 0846+513 and all 
sources within 15° from it was frozen to the value resulting 
from the likelihood analysis over the entire year, while for 
the sources between 15° and 20° from SBS 0846+513 the 
photon index was frozen to the value reported in the 2FGL 
catalogue. If TS < 10 the value of the fluxes were replaced 
by the 2-cr upper limits. The systematic uncertainty in the 
flux is energy dependent: it amounts to 10% at 100 MeV, 
decreasing to 5% at 560 MeV, and increasing to 10% above 
10 GeV. By means of the gtsrcprob tool we estimated that 
the highest energy photon emitted by SBS 0846+513 (with 
probability >80% of being associated with the source) was 
observed on 2 May 2011 at distance of 0.34° from the source 
with an energy of 31 GeV. 

The source was not continuously detected over the en- 
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Figure 2. Integrated flux (E > 100 MeV) light curve of 
SBS 0846+513 obtained from 2011 June 4 to 2011 July 4 with 
3-day or 1-day time bins. Arrows refer to 2-cr upper limits on the 
source flux. Upper limits are computed when TS < 10. 

tire third year and the flux remained below 1 CV 7 ph cm -2 
s _1 , except in the period 2011 June 4-July 4 when an in- 
crease in flux by a factor of ~ 4 was observed. Considering 
the high activity of the source we extracted a spectrum over 
that period, obtaining a photon index F = 1.98 ± 0.05 and 
a flux of (24.4 ± 2.1) xlO -8 photons cnW 2 s _1 . A simi- 
lar spectral evolution has been already observed in other 
fl at spectrum radio quasars (FSRQs) during a bright state 
(Abdo et al.1 l2010dh . A light curve focused on the period 
2011 June 4-, July 4, extracted with the photon index fixed to 
F = 1.98, was built with 3-day and 1-day time bins (Fig. [5]). 
We used a 1-day time bin for the period with higher statis- 
tics. The peak of the emission was observed between June 
30 19:43 UT and July 1 19:43 UT, with a flux of (87 ± 16) 
xlO -8 photons cm -2 s _1 in the 0.1-300 GeV energy range. 
The peak 7 -ray flux is a factor of ~13 higher with respect 
to the average flux estimated over the entire third year of 
observation. The 7 -ray flare is characterized by a sharp in- 
crease in flux of a factor of ~2.5 in 1 day and return to the 
previous flux level in 1 day. 

For taking into consideration the possible influence of the 
choice of the To (the start time of the bin) used for build- 
ing the light curve on the determination of the 7 -ray peak 
we shifted the To of 12 hours back and forth and recalcu- 
lated the daily flux, obtaining (53 ± 15) XlO -8 and (66 ± 
15) xlO -8 photons cm -2 s _1 , respectively. Finally, during 
the month of high activity, replacing the power-law with a 
log parabola or a broken power-law we obtain a TS cur „e 
= 49, indicating a significant curvature. We noted that the 
curvature parameter /3 in the log-parabola spectral model 
increased as the flux increased (see Tabic [TJ. 


3 SWIFT DATA: ANALYSIS AND RESULTS 

The Swift satellite (iGehrels et al.l [20041 ') performed two ob- 
servations of SBS 0846+513 on 2011 August 30 and Septem- 
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Table 1. Unbinned likelihood spectral fit results. 



PL 

LP 


BPL 

Time Period (MJD) 

r 

a 

3 

IT 

r 2 

55412-55777 

2.23 ± 0.05 

1.93 ± 0.12 

0.13 ± 0.05 

2.00 ± 0.09 

2.70 ± 0.18 

55717-55747 

1.98 ± 0.05 

1.23 ± 0.19 

0.30 ± 0.07 

1.51 ± 0.12 

2.81 ± 0.21 


Table 2. Log and fitting results of Swift./X RT observations of SBS 0846+513 using a power-law model with JVh fixed to Galactic 
absorption. “Observed flux. 


Observation 

Net Exposure Time 

Photon index 

Flux 0.3-10 keV“ 

Date 

sec 

r 

XlO -13 erg cm -2 s -1 

2011-08-03 

6738 

1.51 + 0.24 

9.0 ±2.2 

2011-09-15 

5427 

1.43 + 0.28 

7.7 ± 2.9 


ber 15 for 6.8 ksec and 5.5 ksec, respectively. The observa- 
tions were performed with all three on board intruments: the 
X-ray Telescope 1XR.T: [Burrows et al.1 (|2005l ). 0.2-10.0 keV), 
the U ltraviolet Optical Telescope (UVOT; iRoming et al.l 
(120051) , 170-600 r an) a nd the Burst Alert Telescope (BAT; 
iBarthelmv et al J (12005 ). 15-150 keV). 

The hard X-ray flux of this source is below the sensi- 
tivity of the BAT instrument for such short exposures and 
therefore the data from this instrument will not be used. 
Moreover, the source was not present in the Swift BAT 
58-month hard X-ray catalogue (lBaumgartner_et_ah 2010|) 
and t he 54-month Palermo BAT catalogue I Cusumano et all 
120101 ) . 


The XRT data were processed with standard procedures 
(xrtpipeline vO.12.6), filtering, and screening criteria by 
using the Heasoft package (v6.11). The data were collected 
in photon counting mode in both of the observations, and 
only XRT event grades 0-12 were selected. The source count 
rate was low (< 0.5 counts s -1 ); thus pile-up correction was 
not required. Source events were extracted from a circular 
region with a radius of 20 pixels (1 pixel ~ 2.36”), while 
background events were extracted from a circular region 
with radius of 50 pixels away from the source region. An- 
cillary response files were generated with xrtmkarf , and ac- 
count for different extraction regions, vignetting and PSF 
corrections. We used the spectral redistribution matrices 
v013 in the Calibration database maintained by HEASARC. 

Considering the low number of photons collected (< 
200 counts) the spectra are rebinned wit h a minimu m of 
1 count per bin and the Cash statistics dCash I Il979l ) are 
used. We fit the spectrum with an absorbed power-law us- 
ing t he photoelectric absorption model tbabs (IWilms et alJ 
120001 ). with a neutral hydro gen column fixed to its Galac- 
tic value (2.91xl0 20 cm~ 2 : [Kalberla et al.ll2005l ). The low 
photon statistics prevents us from fitting the X-ray data 
with a more complic ated model than a simple power-law. 
We noted that also in lGrupe et alJ (I2OI0I) . in which a large 
sample of NLSls observed by Swift/XRT was studied, the 
absorption column density was fixed to the Galactic value 
and the authors found that in the majority of NLSls the 
spectrum can be fit sufficiently well by an absorbed single 
power-law model with negligible intrinsic absorption. The 
fit results are reported in Table [2] The relatively hard X-ray 


spect rum with respect to the other NLSls ('e.g. lGrupe et all 
l2010l) could be due to the contribution of inverse Compton 
radiation from a relativistic jet, similarly to flat spectrum ra- 
dio quasars. A photon index 1.5-1. 8 was observed in X-rays 
also for PMN J 0948+0022, t he first NLS1 detecte d in 7 rays 
by Eerrm-LAT jAbdo et al.ll2009bl : iFoschini et al.ll201ll ). 

I11 the past SBS 0846+513 was detected in X-rays only 
by ROSAT, with a 0. 1-2.4 keV flux of 2.7xl0~ 13 erg c m~ 2 
s -1 and a photon index T = 1-77+Q gQ ([Yuan et all2008l ). As 
a comparison the flux observed on August 30 by Swift/XRT 
in the same energy range with T fixed to 1.77 is 3.7 ± 0.6 
xlO -13 erg cm -2 s _1 , ~40% higher than the ROSAT obser- 
vations, adding evidence of the variability character of the 
source. 

In 2011 August the UVOT instrument took 5 frames in 
the wl band only, while four sequences with all v, b, u, wl, 
m2, and w2 filters were acquired in 2011 September. These 
data were processed with v6.10 of the Heasoft package and 
the CALDB release dated 2011 August 12. Source counts 
were extracted from a circular region of 5 " radius centred 
on the source, while background counts were derived from a 
circular region of 10” radius in the source neighbourhood. By 
fitting the source spectrum of September 15 with a power- 
law, we calculated the effective wavelengths, count rate to 
flux conversion factors (CFa), and Galactic extinctions for 
the UVOT bands accord ing to the procedure explained in 
iRaiteri et al.l d2010l . [201 ll) . The results are shown in Table O 
To obtain the UVOT de-reddened fluxes reported in Table [4] 
we multiplied the count rates for the CFa and corrected for 
the corresponding Galactic extinction values Aa- 

We note that a Swift/UVOT observation performed on 
2011 September 15 found SBS 0846+513 about 0.4 magni- 
tude brighter in V-band with respect to the HST observation 
in January 1992. 


4 RADIO DATA: ANALYSIS AND RESULTS 

New and archival radio data for SBS 0846+513 were col- 
lected from Medicina, OVRO, the Effelsberg 100-m, VLA 
and VLB A. In the following subsections we present the ob- 
servations performed by these facilities. 
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Table 3. Results of the UVOT calibration procedure: effec- 
tive wavelengths A e fj, count rate to flux co nversion factors CF a , 
and Galactic extinction calculated from the lCardelli et ahl lil989h 
laws. 


Filter 

-Vff 

A 

CF A 

10 -16 erg cm -2 s -1 A -1 

Aa 

mag 

V 

5430 

2.60 

0.09 

b 

4360 

1.47 

0.11 

u 

3476 

1.65 

0.14 

uvw 1 

2621 

4.40 

0.19 

uvm2 

2257 

8.36 

0.24 

uvw 2 

2087 

5.97 

0.23 


Table 4. Results of the analysis of UVOT data for SBS 0846+513. 


MJD 

Filter 

Magnitude 

Magnitude.error 

55803.644 

UVW 1 

19.57 

0.07 

55819.497 

V 

19.10 

0.30 

55819.489 

b 

19.67 

0.19 

55819.488 

u 

19.22 

0.17 

55819.497 

uvwl 

19.41 

0.15 

55819.501 

uvm2 

19.55 

0.15 

55819.493 

uvw2 

19.55 

0.10 


4.1 Medicina 

We observed SBS 0846+513 with the Medicina radio tele- 
scope five times between 2011 August and November. We 
used the new Enhanced Single-dish Control System (ESCS) 
acquisition system, which provides enhanced sensitivity and 
supports observations with the cross scan technique. All ob- 
servations were performed at both 5 and 8.4 GHz; the typical 
on source time is 1.5 minutes and the flux density was cali- 
brated with respect to 3C 286. Since the signal to noise ratio 
in each scan across the source was low (typically ~ 3), we 
performed a stacking analysis of the scans, which allowed 
us to significantly improve the signal to noise ratio and the 
accuracy of the measurement. We list the final values of the 
5 and 8.4 GHz flux density in Table [S] 

4.2 OVRO 

As part of an ongoing blazar monitoring program, the Owens 
Valley Radio Observatory (OVRO) 40 m radio telescope has 
observed SBS 0846+513 at 15 GHz regularly since the end 
of 2007 llRichards et al.|[201ll) . This monitoring program in- 
cludes over 1500 known and likely 7 -ray loud blazars above 

Table 5. Results of the Medicina 32-m radio observations. 


Obs. date 

s 5 

(mjy) 

Ss.4, 

(mjy) 

2011-08-10 

170 ± 15 

235 ± 20 

2011-09-08 

210 + 20 

185 ± 20 

2011-09-22 

190 ± 20 

180 ± 20 

2011-10-13 

190 ± 20 

280 ± 30 

2011-11-16 

250 ± 30 

275 ± 25 


2009 201 0 2011 2012 



Time [ MJD] 

Figure 3. 15 GHz radio light curve for the period 2008 August 
4-2011 December 12 from the OVRO telescope. The downward 
arrow indicates the time of the peak of the 7 -ray activity observed 
by Fermi- LAT. 

declination —20°. The sources in this program are observed 
in total intensity twice per week with a 4 mjy (minimum) 
and 3% (typical) uncertainty. Observations are performed 
with a dual-beam (each 2.5 arcmin FWHM) Dicke-switched 
system using cold sky in the off-source beam as the refer- 
ence. Additionally, the source is switched between beams 
to reduce atmospheric variations. The absolute flux density 
scale is calibrated using obser vations o f 3 C 286, adopting 
the flux density (3.44 Jy) from iBaars et al.l (ll977l) . This re- 
sults in about a 5% absolute scale uncertainty, which is not 
reflected in the plotted errors. 

4.3 EfFelsberg 100-m 

The centimeter spectrum of SBS 0846+513 was observed 
with the Effelsberg 100-m telescope on 2011 April 30 (MJD 
55681.8) within the framework of a Fermi - related mon- 
itoring program of 7 -ray blazars (F-GAMMA program; 
IFuhrmann et aI.M2007f ). The measurements were conducted 
with the secondary focus heterodyne receivers at 2.64, 8.35, 
14.60, and 32.00 GHz. The observations were performed 
quasi-simultaneously with cross-scans, that is, slewing over 
the source position, in azimuth and elevation directions, with 
adaptive numbers of sub-scans for reaching the desired sensi- 
tivity (for details, see lFuhrmann et alJboQsl i lAngelakis et all 
120081) . Pointing offset correction, gain correction, atmo- 
spheric opacity correction, and sensitivity correction have 
been applied to the data. 

4.4 VLA and VLBA data 

To implement the information on the flux density variability 
available from the single-dish observations and to study the 
source structure on kpc and pc scales, we analysed archival 
VLA and VLBA data at different frequencies obtained be- 
tween 1986 and 2011. To study the proper motion we also 
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Table 6. Log of the archival VLA observations and flux density. 


Freq 

GHz 

Date 

Code 

Obs.time 

min 

Beam 

arcsec x arcsec 

Flux 

mjy 

1.4 

1986-04-10 

AV 127 

9.16 

1.40x1.12 

194+6 

4.8 

1986-04-10 

AV 127 

9.16 

0.44x0.35 

286+9 

4.8 

1996-01-05 

BW 021 

1.66 

1.39x1.19 

332+10 

4.8 

1996-12-30 

baoi8 

3.33 

0.43x0.36 

363+11 

4.8 

2009-05-09 

CALSUR 

0.66 

2.12x1.18 

196+6 

8.4 

1995-08-13 

AM484 

0.5 

0.24x0.22 

356+11 

8.4 

1995-09-02 

AM484 

0.5 

0.32x0.23 

327+10 


made use of 15-GHz VLBA data from the MOJAVE pro- 
grammed 

Usually SBS 0846+513 was observed with the VLA as a 
phase calibrator and the on-source observing time is gener- 
ally quite short. Logs of VLA and VLBA observations are 
reported in Tables [6] and [7] respectively. 

The data reduction of both VLA and VLBA data was per- 
formed following the standard procedures implemented in 
the National Radio Astronomy Observatory (NRAO) AIPS 
package. In the case of VLA data, the accuracy of the am- 
plitude calibration was checked by means of either 3C48 
or 3C 286, and it resulted to be within 3%. For the VLBA 
datasets, a priori amplitude calibration was derived using 
measurements of the system temperature and the antenna 
gains. Uncertainties of the amplitude calibrations were esti- 
mated to be within 5-10%. For MOJAV E data we imported 
the fully calibrated uv datasets dLister et alj|2009h . 

The final VLA and VLBA images were produced after a 
number of phase-only self-calibration iterations. The source 
parameters, like flux density and angular size, were derived 
on the radio image plane by means of the AIPS task JMFIT 
which performs a Gaussian fit to the source and its subcom- 
ponents (Table [8}. For the extended component visible at 15 
GHz, we determine the flux density by means of TVSTAT, 
which performs an aperture integration on a selected region 
on the image plane. 


5 RADIO VARIABILITY 

Multi-epoch studies of the radio emission of SBS 0846+513 
show significant flux density variability. In Fig. [3] we report 
the OVRO observations obtained since the launch of Fermi. 
From the analysis of the OVRO light curve at 15 GHz it 
is clear that high and low activity states are interspersed. 
Interestingly, no 7-ray emission was detected during the ma- 
jor radio flaring episodes that occurred before 2011. A max- 
imum flux density variation of a factor of ~2.0 was observed 
over 3 months during the period 2010 August-2011 August, 
while in the period 2008 August-2010 August a variation of 
a factor of ~2.8 was observed over a similar time interval. 
Flux density variability is detected also at 4.8 GHz, where 
the flux density varies about 70% during the time inter- 
val spanned by the VLA data (1986 to 2009). However, the 
poor time coverage of the VLA data does not allow us to 

4 The MOJAVE data archive is maintained at 
http: / /www. physics. purdue.edu/MOJAVE 


2011.6 201'. 65 2011.7 2011.75 2011.8 



Figure 4. Radio data of SBS 0846+513 collected at 5.0 GHz 
(triangles) and 8.4 GHz (squares) by Medicina, and 15 GHz (open 
circles) by OVRO during the period 2011 July 28 - November 19 
(MJD 55770-55884). 


fully characterize the flux density behaviour on such a long 
time interval. 

The detection of the 7-ray flare in 2011 June, triggered a 
monitoring campaign with the Medicina telescope. In Fig. [4] 
we report the Medicina flux density at 5 ( triangles ) and 8.4 
GHz ( squares ) in addition to the OVRO 15-GHz data ( open 
circles) collected between 2011 July 28 and November 19. 
From August 20 (MJD 55793.628), i.e. about one month af- 
ter the 7-ray flare, the flux density at 15 GHz started to 
increase reaching the maximum of about 280 mjy around 
October 18, and then it started to decrease. At lower fre- 
quencies the flux density variation is observed after a longer 
time, firstly at 8.4 GHz (October 13), and then at 5 GHz 
(November 16). This time delay may be explained by opacity 
effects, which are more severe at longer wavelengths. Dur- 
ing this time interval also the shape of the radio spectrum 
changes, as is clearly visible in Fig. 0 (left panel). We must 
note that Medicina and OVRO observations are not strictly 
simultaneous and to build the radio spectrum we considered 
the OVRO observation that is the closest in time to the 
Medicina observing epoch. For comparison, in Fig. [5] (right 
panel) we show the simultaneous radio spectrum between 
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Figure 5. Left panel : Radio spectra of SBS 0846+513. Observations at 5 and 8.4 GHz are from the Medicina single-dish telescope, those 
at 15 GHz are from OVRO. Triangles, squares, and crosses refer to observations performed around 2011 August 10, October 13, and 
November 16, respectively. Right panel: Radio spectrum of SBS 0846+513 obtained by the Effelsberg single-dish on 2011 April 30 from 
2.64 to 32 GHz. 


VLA 8.4 GHz 



Figure 6. VLA image at 8.4 GHz of SBS 0846+513. On the im- 
age we provide the peak flux density, in mjy/beam, and the first 
contour intensity (f.c., in mjy/beam) that corresponds to three 
times the noise measured on the image plane. Contour levels in- 
crease by a factor of 2. The beam is plotted on the bottom left 
corner of the image. 

2.64 and 32 GHz obtained by Effelsberg on 2011 April 30, 
i.e. before the high activity state detected in 7 rays. 


6 RADIO MORPHOLOGY 

The resolution of the VLA is not adequate to resolve the 
radio structure of SBS 0846+513 (Fig. [6j even when the 
array is in its most extended configuration (full width 


half maximum ^0.3 arcsec). When imaged with the high 
spatial resolution of the VLBA the source is resolved in two 
components with a core-jet structure (Figs. [T]_and El), as 
also pointed out in a previous work by I Taylor et al.l J2005h . 
The core (component W) is still unresolved with an upper 
limit of 0.3 mas, while the jet structure (component E) is 
3.5 x 1.6 mas. The flux density ratio between components 
E and W is ~19 and 14 at 5 and 8.4 GHz respectively. 
Their separation is about 3.5 mas, corresponding to ~23 pc 
projected distance, given the source’s redshift. At 15 GHz 
component E shows an extended low-surface brightness 
structure with a steep spectrum. On the other hand 
component W is resolved into two compact components 
(labelled W1 and W2 in Fig. |9j). To investigate a possible 
proper motion of the jet, we compared the separation 
between Wl, considered the core region, and W2, assumed 
to be a knot in the jet, at the three observing epochs. To 
this purpose, in addition to the data analysis on the image 
plane, we model-fitted the visibility data using gaussian 
components of the three-epoch MOJAVE data by means 
of the model- fitting option in difmap. This approach is 
preferable in case we want to derive small variations in 
the source structure and little changes in the position of 
the source components. From this comparison we found 
that Wl and W2 are separating with an apparent velocity 
of (0.38±0.07) mas/yr (Fig. 1101) . which corresponds to 
(8.2±1.5)c. This apparent superluminal velocity suggests 
the presence of boosting effect. However, the availability of 
only three observing epochs spanning a short time baseline 
without frequent time sampling implies large uncertainties 
on the estimated values. 

When we compare the total flux density derived from 
VLBA data with that from the VLA, we find that the 
VLBA can recover at most 85% of the VLA flux density, 
even when the observations were performed simultaneously 
(i.e. 1996 January 6). This may be related to the slightly 
different observing frequencies, i.e. 4.8 GHz at the VLA, 
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Figure 7. VLBA image at 5.0 GHz of SBS 0846+513. On the 
image we provide the peak flux density, in mJy/beam, and the 
first contour intensity (f.c., in mJy/beam) that corresponds to 
three times the noise measured on the image plane. Contour levels 
increase by a factor of 2. The beam is plotted on the bottom left 
corner of the image. 


Figure 8. VLBA image at 8.4 GHz ( bottom ) of SBS 0846+513. 
On the image we provide the peak flux density, in mJy/beam, and 
the first contour intensity (f.c., in mJy/beam) that corresponds 
to three times the noise measured on the image plane. Contour 
levels increase by a factor of 2. The beam is plotted on the bottom 
left corner of the image. 


Table 8. Flux density in mjy of the components of SBS 0846+513 
from VLBA data. Epochs of the observations are reported in Ta- 
ble [7] 



E 

W 

W1 

W2 

S5GHz,epl 

143=2 

2693=19 

- 

- 

S5GHz,ep2 

17±2 

2653=18 

- 

- 

38.4GHz 

14±2 

1963=14 

- 

- 

Sl5.3GHz,epl 

3.43=0.4 

1923=15 

1403=10 

523=4 

Sl5.3GHz,ep2 

3.1±0.3 

2403=24 

1753=17 

653=6 

SlS.3GHz,ep3 

3.6±0.4 

2093=20 

1513=15 

583=5 


whereas it is 5.0 GHz at the VLBA. However, in this case 
the spectral index would have to be extremely steep (a r ^4; 
f u oc u~ ar ) to obtain such a difference in the flux density, 
making this explanation unlikely. Another possibility is 
related to observational limitations due to the lack of short 
spacings of the VLBA array, implying that only structures 
smaller than ~40 and 30 mas at 5 and 8.4 GHz, can be 
detected. This suggests that the missing flux on the parsec 
scale may be related to extended, low-surface brightness 
features like a jet component resolved out by the VLBA 
array. 


VLBA 15.3 GHz 
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7 SED MODELING 

The 7 -ray detection of the radio loud NLSls has given us 
the possibility of studying the characteristics of this class of 
objects by modeling their broad band spectra. 

We have built a non- simultaneous SED, all of which is 


Figure 9. VLBA image at 15.3 GHz of SBS 0846+513. Data 
are from the MOJAVE programme. On the image we provide the 
peak flux density, in mJy/beam, and the first contour intensity 
(f.c., in mJy/beam) that corresponds to three times the noise 
measured on the image plane. Contour levels increase by a factor 
of 2. The beam is plotted on the bottom left corner of the image. 
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Table 7. Log of the VLBA observations analysed in this paper and total flux density. a : data from the MOJAVE programme. 


Freq 

GHz 

Date 

Code 

Beam 
mas x mas 

Obs time 
min 

Flux 

mjy 

5.0 

1995-03-15 

BW 015 

3.99x1.34 

50 

284±20 

5.0 

1996-01-06 

BW 021 

4.13x1.29 

50 

281±20 

8.4 

2011-03-12 

BC196J 

2.00x0.97 

10 

210±16 

15.3 

2011-05-26 

RL149DT" 

0.77x0.59 

30 

244T17 

15.3 

2011-07-15 

BL149DM a 

0.76x0.60 

30 

213±15 

15.3 

2012-01-02 

BL178AE a 

0.74x0.64 

36 

231±21 



Year 


Figure 10. Changes in separation with time between components 
W1 and W2. The solid line represents the regression fit to the 15- 
GHz MOJAVE data. 


based on non-flaring data. This ought to adequately repre- 
sent the emission from this object in a low state. The LAT 
spectrum was built with data centred on 2010 October 4 to 
2011 June 4 (MJD 55473-55716). In addition we included 
in the SED the Effelsberg radio data collected on 2011 April 
30 and the Swift (UVOT and XRT) data collected on 2011 
September 15, thus well before and after the 7-ray flare. The 
data from 2MASS All-sky PSC source catalogue and WISE 
preliminary release Source Catalogue provided information 
about the IR part of the spectrum. The flare centred on ~ 
MJD 55750 had a variability timescale of ~ 1 day, which 
constrains the size of the emitting region during the flare. 

We modeled the SED assuming emission from a rel- 
ativistic jet with mechanisms of synchrotron, synchrotron 
self-Compton (SSC), and Compton-scattering of a dust 
torus external to the jet (EC-dust ) . The des cription of the 
model can be found in I Finke et al.l d2008l J and lDermer et abl 
(2009). Jet powers were calculated assuming a two-sided jet. 
Although the flaring timescale does not constrain the emit- 
ting region during the low state, we have produced a model 
fit which is roughly consistent with this timescale (about 1 
day). The synchrotron component can adequately explain 
the IR. through UV points, although some of the 2MASS 
points are significantly higher than the rest, and are not 
well- fit. This is probably due to contamination by IR emis- 
sion from the host galaxy. We note that this IR excess could 
be consistent with the typical starlight expected from an el- 
liptical galaxy. The radio spectrum is flat in F v (Fig.O and 
probably the result of a superposition of several jet compo- 


nents (e.g., iKonigl Ill98ll '). so for our purposes these points 
are considered upper limits. The synchrotron component we 
use is self-absorbed below ~ 10 12 Hz. An attempt was made 
to fit the X-ray through 7 ray data with an SSC component, 
but this was not found to be possible. Instead, an external 
component was required. Correlations of 7-ray and optical 
flares with radio light curves and rotations of optical polar- 
ization angles in low-synchrotron-peaked blazars seem to in- 
dicate the 7-ray /optical emitting region is outside the broad 
line region, where the dust torus is the likely seed photon 
source fe.g.. iMarscher et al.ll2C)foh . So we chose as our seed 
photon source a dust torus, which was simulated as a one- 
dimensional ring with radius Rdust. aligned orthogonal to the 
jet, emitting as a blackbody with temperature Tdust and lu- 
minosity Ldust- Observational evidence for a dust torus in 
NLSls as well as in the othe r Seyfert galaxies was recently 
reported in lMor et al.l (l2012h . The fit can be seen in Fig. [TT] 
and the parameters can be found i n Ta ble [HI _A descr iption 
of the parameters can be found in iDermer et al.l (|2009| J . The 
dust parameters were chosen so that Rdust is rou ghly con - 
sistent with the sublimation radius dNenkova et alj|2008h . 
assuming the torus luminosity is about 1/10 of the disc lumi- 
nosity. If the emitting region takes up the entire width of the 
jet, and the jet is assumed to be conical, its half-opening an- 
gle would be about 1°, approximately consist ent wi th those 
meas ured from VLBI observations of blazars (Ijorstad et al.l 
l2005ll . This size scale is roughly the same as for the 15 GHz 
core (Fig. [9j, and so if this model fit is correct, the 7-ray 
emitting region is probably located near this region, which 
is presumably a synchrotron self-absorption photosphere. 

The electron distribution used, a broken power-law with 
index pi = 2.2 below the break at Ark and P2 = 3.2 for 
7 < A: is consistent with particles injected with index 
2.2, and emission taking place in the slow cooling regime 
(e.g., iBottcher fe Dermerll20oT That is, particles are in- 
jected between Amin and Amax, with a cooling break at 

3 m e c 2 

'Ibrk = ~ 7 

^CO’T’tesc'U'tot 

where u tot. is the total energy density in the blob frame, 
which in the case of our model fit is dominated by the ex- 
ternal energy density. The escape time t esc ~ R'b/ C so that 
for our model fit 7 ' brk ~ 300. Also note that in this model fit 
the magnetic field and electrons are nearly in e quipartitio n. 

Using the SDSS Data Release (DR) 1, Izhou et ahl 
d2005lJ estimated the black hole mass from several methods, 
the H/3 broad line, the [OIII]A5007 narrow line, and the host 
galaxy’s bulge luminosity, as ~ 8.2 x 10 6 Mq,~5.2x10 7 Mq, 
and ~ 4.3 x 10 ' Mq, respectively. Large uncertainties (0.4- 
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Figure 11. Spectral energy distribution data (circles and 
squares) and model fit (solid curve) of SBS 0846+513 with the 
model components shown as dashed curves. The data points were 
collected by Effelsberg (2011 April 30), Swift (UVOT and XRT; 
2011 September 15), and Fermz-LAT (2010 October 4-2011 June 
4), together with archival data from 2MASS and WISE. The LAT 
spectrum during the flaring period (2011 June 4- July 4) is shown 
as the diamonds. 


0.7 d e x) are asso c iated to these techniques (jVerstergaardl 
120041 ). S hen et all (| 20 111 ) estimated the BH mass based on 
the H,9 broad line and Mgll line in the SDSS DR7 spectrum 
as log(A/sir /Mq) = 7.99 ± 0.12 and log(Msir /Mq) = 7.79 
± 0.16, respectively. Note that the broadband SED shows no 
evidence for a blue bump, so the optical spectrum probably 
has considerable contamina tion Jrom jet emission not taken 
into account in IShen et al.l (1201 lh . and s o these mass esti - 
mates should be taken with caution. If the lZhou et al.l (l2005l l 
H/l estimate is correct, the total jet power will exceed the 
Eddington luminosity for this source ( Lsdd = 1.3 x 10 45 
erg s _ 1 Mbh / ( 10 7 M@) ). Thus the BH mass is probably 
more on the high side of these estimates, if the jet power is 
Eddington- limited. 

The Compton dominance for SBS 0846+513, i.e., the 
ratio of the peak luminosities of the Compton and syn- 
chrotron components, is « 7, which is a rather standard 
value for FSRQs, although quite a high value for BL Lacs 
(e.g., ICiommi et alJ lioill : iFinke I [201+ . Its SED in many 
ways resembles that of a FSRQ, with an X-ray spectral in- 
dex ~ 1.5, consis tent with values for FSRQs in the BAT 
(lAiello et a.l .11200?!) and BeppoSAX (iDonato et alJ2005h cat- 
alogues. As we discuss in section [5] its LAT flux and spectral 
index are also consistent with values for FSRQs. 

The flaring LAT spectrum is also shown in Fig. 1111 Un- 
fortunately, there are no multiwavelength data simultaneous 
to this flare. For FSRQs, 7-ray flares are often associated 
with optical flares, although there are occasions when they 
are not. Both types of flares are seen, for ex ample, from the 
FSRQ PKS 1510—089 dMarscher et al.l[2010l ) so we can only 
guess at the behavior of SBS 0846+513 during this flare. 
However, note that during the flare, the object reaches a 
luminosity of I/ 7 ,i so ~ 1.0 x 10 48 erg s _1 , making the radia- 
tive power L 7tra d ~ L 7 , iso /r 2 ~ 4.4x 10 45 erg s+ assuming 
r = 15. This is about 1/2 of the Eddington luminosity for a 
10 7 Mq black hole. 


8 DISCUSSION AND CONCLUDING 

REMARKS 

After the 4 o bjects detect ed by Fermi - LAT in the first year 
of operation dAbdo et alj|2009cl h SBS 0846+513 is a new 
NLS1 detecte d by Fermi - LAT during high 7-ray activity 
in 2011 June dDonato fe Perki~nsll2011 h The power released 
by this object during the flaring activity was a strong in- 
dication of the presence of a relativistic jet as powerful as 
those in blazars, supported by the apparent superluminal 
velocity of the jet derived by tracking the position of a jet 
component in 2011-2012 VLBA data (see Sec. [6]). Before 
the 7-ray flaring episode, the simultaneous multifrequency 
observations performed by Effelsberg showed a flat radio 
spectrum up to 32 GHz. After the flare, the spectral shape 
changed, becoming convex. The spectral variability was also 
accompanied by variations in the radio flux density, which 
were originally detected at the higher frequencies, later mov- 
ing to lower frequencies, likely due to opacity effects. These 
spectral and variability properties indicate blazar-like be- 
haviour ! _nduch_Jins_a^ady_been observed in other 7-ray 
NLSls IIFuhrmann et al1l201ll b 

iGhisellini et aid d2008l . l201ll L investigating the 7-ray 
properties of blazars detected by Fermi - LAT in the first year 
of observation, suggested a transition between BL Lac ob- 
jects and FSRQs that can be justified mainly by the different 
accretion regimes: highly sub-Eddington in the former case, 
near-Eddington in the latter. In this context the high accre- 
tion rate of SBS 0846+513, and more generally of the 7-ray 
radio- loud NLSls, is another indication of its similarity with 
FSRQs. A comparison of the BL L acs and FS RQ s fro m the 
First LAT AGN Catalog (1LAC; lAbdo et al.ll 20101 ) with 
“misaligned AGN” detected by Fermi (MAGN; non-blazar 
AGN with jets pointed away from the observer) in the F 7 - 
L 7 plane has shown that MAGN and blazars occupy differ- 
ent regions of the plane, with only two high redshift FRII 
galaxies, 3C 207 and 3C 380, which lie among the FSRQs 
IIAbdo et alJl2010bh . This is in agreement with the idea that 
MAGN are less beamed than blazars. In this context it is 
interesting to consider also the 7-ray NLS ls. For a direct 
comparison with the results shown in lAbdo et al.l (|2010bl ~) 
we calculated the flux and photon index over the third year 
of Fermi operation in the 0.1-10 GeV energy band, resulting 
in F 7 = 2.19 ± 0.06 and Flux (0.1-10 GeV) = (6.6+0.6) 
x 10~ 8 erg cm -2 s _1 . The corresponding observed isotropic 
7-ray luminosity in the 0.1-10 GeV is 3.6xl0 46 erg s _1 . We 
plotted these values of SBS 0846+513 in the F 7 - L 7 plane 
toge ther with th e FSRQs, BL Lacs, and misaligned AGN 
from lAbdo et alJ (l2010bl L As can be seen SBS 0846+513 lies 
in the blazar-region, in particular in the transition region be- 
tween the distribution of BL Lacs and FSRQs (Fig. 1121) . We 
note that also the other 7-ray NLS1 observed in flaring activ- 
ity, PMN J0948+0022 (F 7 = 2.26 ± 0.08 and 0.1-10 GeV 
luminosity of 9.6 xl0 4 6 erg s~ 4 over the first 24-month of 
Fermi operation: lGrandI1l201ll L occupies the same blazar- 
region in that plane. This should reflect a similar viewing 
angle with respect to the jet axis and beaming factors for 
the 7-ray emission between blazars and the two 7-ray NLSls 
SBS 0846+513 and PMN J0948+0022. In the same way the 
spectral evolution during the flaring activity in June 2011 
observed in 7 rays from SBS 0846+513 is a common be- 
haviour in bright FSRQs and low-synchrotron-peaked BL 
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Table 9. Model parameters for the SED shown in Fig, fill 


Redshift 

z 

0.5835 

Bulk Lorentz Factor 

r 

15 

Doppler Factor 


15 

Magnetic Field 

B 

1.0 G 

Variability Timescale 

tv 

lx 10 s s 

Comoving radius of blob 

K 

2.8xl0 16 cm 

Jet Height 

r 

1.6 x 10 18 cm 

Low-Energy Electron Spectral Index 

pi 

2.2 

High-Energy Electron Spectral Index 

P2 

3.2 

Minimum Electron Lorentz Factor 

y ■ 

1 min 

5.0 

Break Electron Lorentz Factor 

Ihrlc 

3.0 x 10 2 

Maximum Electron Lorentz Factor 


9.0 x 10 3 

Dust Torus luminosity 

■t'dust 

4.0 X 10 44 erg s _1 

Dust Torus temperature 

Tdust 

1.5 x 10 3 K 

Dust Torus radius 

R-dust 

2.0 X 10 18 cm 

Jet Power in Magnetic Field 

Pj,B 

1.4 X 10 45 erg s _1 

Jet Power in Electrons 

P- 

1 j,par 

4.3 X 10 44 erg s _1 


Lacs detected by Fermi (|Abdo et al.ll2010dh . with a change 
in photon index < 0.2-0. 3 and an increasing spectral curva- 
ture. 

One of the key questions is the maximum power re- 
leased by the jets of NLSls. During the 7-ray flaring ac- 
tivity observed in 2011 June-July SBS 0846+513 reached 
an observed isotropic 7-ray luminosity (0.1-300 GeV) of 
l.OxlO 48 erg s 1 on daily timescales, compar able to that of 
luminous FSRQs. After PMN J0948+0022 dFoschini et alJ 
B), this is the second NLS1 observed to generate such a 
high power. This could be an indication that all the radio- 
loud NLSls are able to host relativistic jets as powerful 
as those in blazars, despite the lower BH mass; alterna- 
tively some NLSls could have peculiar characteristics al- 
lowing the development of these relativistic jets. The mech- 
anism at work for producing a relativistic jet is not clear, and 
the physical parameters that drive jet formation is still un- 
der debate. One fundamental parameter could be the black 
hole mass, with only large masses allowing for the efficient 
formation of a rel ativistic jet. It was fo r example noted 
by some authors (iMcLure fc Jarvij 120041 : iLiu et all 120061 : 
ISikora et al.ll2007h that quasars with Mbh > 10 8 Mq reach 
radio loudness 3 orders of magnitude greater than quasars 
with A/bh < 3 x 10 7 Mq. The large radio loudness of SBS 
0846+513 could chall enge this id ea if the black hole masses 
estimated by IZhou et al.l J2005I) are confirme d. According 
to the “modified spin paradigm” discussed in ISikora et al.l 
(120071') . another fundamental parameter for the efficiency of 
a relativistic jet production should be the BH spin, with 
SMBHs in elliptical galaxies having on average much larger 
spins than SMBHs in spiral galaxies. This is due to the fact 
that spiral galaxies are characterized by multiple accretion 
events with random orientation of angular momentum vec- 
tors and small increments of mass, while elliptical galaxies 
underwent at least one major merger with large matter ac- 
cretion triggering an efficient spin up of the SMBH. The 
accretion rate (thus the mass) and the spin of the black hole 
may be related to the host galaxy, leading to the hypothe- 
sis t hat rel ativi stic je ts can develop only in elliptical galaxy 
(e.g. lMarschei~ll2009l'l . 

I11 this context the discovery of relativistic jets in a 


class of AGN usually hosted by spiral galaxies, such as 
the NLSls, was a great surprise. Unfortunately only very 
sparse observations of the host galaxy of radio- loud NLSls 
are available at this time, and the sample of objects studied 
by iDeo et ahl d2006l ) and IZhou et all (I2OO6I) had z < 0.03 
and z < 0.1, respectively, including both radio-quiet and 
radio-loud objects. We note that BH masses of radio-loud 
NLSls are generally larger with respe ct to the entire sampl e 
of NLSls (A/bh ~ (2 — 10) x 10 7 Af Q iKomossa et alj|200q ). 
even if still small when compared to radio-loud quasars. The 
larger BH masses of radio-loud NLSls with respect to radio- 
quiet NLSls could be related to prolonged accretion episodes 
that can spin-up the BHs. The small fraction of radio-loud 
NLSls with respect to radio- loud quasars could be an indi- 
cation that in the former the high accretion usually does not 
last s ufficiently long to significantly spin-up the BH (ISikora I 
l2009h . 

Of the 7-ray emitting NLSls, host galaxy imaging is 
available only for 1H 0323+342, with HST and the Nordic 
Optical Telescope. These observations reveal a one-armed 
galaxy morphology or a circumnuclear ring, sug gesting two 
possi bilities: the spiral arms of the host galaxy (Zhou_e£_al. 
2007) or the residual of a merging galaxy (lAnton et al . 
200a ). No significant resolve d structures have been observed 
by HST for SBS 0846+513 dMaoz et al.lll993l ). and no high- 
resolution observations are available for the remaining 7-ray 
NLSls. Thus the possibility that the development of rela- 
tivistic jets in these objects could be due to strong merger 
activity, is not ruled out. Further high-resolution observa- 
tions of the host galaxies of 7-ray NLSls will be fundamental 
to obtain important insights into relativistic jet formation 
and development. 

To conclude, SBS 0846+513 shows all the characteris- 
tics of the blazar phenomenon. The extreme power released 
by SBS 0846+513 during the high 7-ray activity in 2011 
July confirms that, as with PMN J0948+0022, NLSls can 
host relativistic jets as powerful as blazars. Radio and 7- 
ray data collected for SBS 0846+513 suggest spectral and 
variability properties similar to blazars, and the modeling of 
the average SED gives similar results to those of blazars, in- 
cluding similar Lorentz factors. This could be an indication 
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Log 1^ [100 MeV-10 GeV] (erg s“‘) 


Figure 12. The photon index T 7 of FRI radio galaxies (red cir- 
cles), FRII radio galaxies (green squares), BL Lacs (open blue cir- 
cles) and FSRQs (open black squares) are plotted together with 
the NLS1 SBS 0846+513 (magenta point including a black trian- 
gle) as a function of the obser ved isotropic 7- ray l uminosity (100 
MeV-10 GeV). Adapted from lAbdo et alJ (1201011 ). 

that these 7 -ray NLSls are low mass (and possibly younger) 
version of blazars. The detection of new radio-loud NLSls 
in 7 rays by Ferrm-LAT will be important for extending the 
source sample, and better characterizing this new class of 
7 -ray emitting AGN. Equally important will be to perform 
further multifrequency observations of the 7 -ray emitting 
NLSls already detected by Fermi and investigate in detail 
their characteristics over the entire electromagnetic spec- 
trum, to help understand their nature. 
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